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Abstract 

A density operator of a bipartite quantum system is called robustly 
separable if it has a neighborhood of separable operators. Given a bi- 
partite density matrix, its property to be robustly separable is reduced, 
using the continuous ensemble method, to a finite number of numerical 
equations. The solution of this system exists for any robustly separa- 
ble density operator and provides its representation by a continuous 
mixture of pure product states. 

1 Introduction 

Entanglement is treated as a crucial resource for quantum computation. It 
plays a central role in quantum communication and quantum computation. 
A considerable effort is being put into quantifying quantum entanglement. 
Usually the efforts are focused on quantifying entanglement itself, that is, 
describing the impossibility to prepare a state by means of LOCC (local 
operations and classical communications). One may, although, go another 
way around and try to quantify separability rather than entanglement: this 
turned out to be applicable for building combinatorial entanglement patterns 
for multipartite quantum systems |2j. Briefly, the contents of this paper is 
the following. 

I dwell on the case of bipartite quantum systems. A state of such system 
is called separable if it can be prepared by LOCC, and ROBUSTLY sepa- 
rable if it has a neighborhood of separable operators in the space of self 
adjoint operators. In terms of density matrices that means that p, its density 
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matrix, can be represented as a mixture of pure product states. I suggest 
to replace finite sums of projectors by continuous distributions on the set of 
unit vectors. The case of single particle system is considered in section |2 
as an example, the explicit form for the equations for the parameters of the 
mixture is obtained. 

In section |21 the bipartite case is considered. The density operators are 
represented by distribution on the Cartesian product of unit spheres in sub- 
systems' spaces. Given a density operator p, we consider it as an element 
of the space £ of all self adjoint operators in the product space 23 = Ti ® Ti. 
Then the robust separability of p is equivalent (see also 4 ) to the solvability 
of the following vector equation in £: 

VK = p 

for the trace functional K on £ whose form is obtained in section|21 When we 
fix a product basis in 53, this equation becomes a system of n transcendent 
equations with respect to n 4 variables. 

2 Smeared spectral decomposition and optimal en- 
sembles 

The set of all self-adjoint operators in Ti = C n has a natural structure of 
a real space M 2n , in which the set of all density matrices is a hypersurface, 
which is the zero surface of the affine functional TrX — 1. 

Generalizing the fact that any convex combination of density operators 
is again a density operator, we represent density operators as probability 
distributions on the unit sphere in the state space Ti of the system. 

Technical remark. Pure states form a projective space rather than the 
unit sphere in Ti. On the other hand, one may integrate over any proba- 
bilistic space. For technical reasons I prefer to represent ensembles of pure 
states by measures on unit vectors m. Ti. I use the Umegaki measure on 
CB n — the uniform measure with respect to the action of U(n) normalized 
so that 

dS n = 1 (1) 

4>eCB n 

Similarly, for bipartite case the integration will be carried out over the Carte- 
sian product of unit spheres in appropriate state spaces. 
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Now pass to a more detailed account of this issue beginning with the 
case of a single quantum system. Let TL = C n be a n-dimensional Hermitian 
space, let p be a density matrix in TL. We would like to represent the state 
whose density operator is p by an ensemble of pure states. Let this ensemble 
be continuous with the probability density expressed by a function p(cj)) 
where <j) ranges over all unit vectors in TL (see the technical remark above). 
The density operator of a continuous ensemble associated with the measure 
p((j)) on the set CB n of unit vectors in TL is calculated as the following 
(matrix) integral 



P 



j p(0) 1^(01 dS n (2) 



where |0)(0| is the projector onto the vector (</>|. Effectively the operator 
integral p @ can be calculated by its matrix elements. In any fixed basis 
{|ej}} in TL, each matrix element pij = (ej|p|e,-) is the following numerical 
integral: 



Pij = (ei\p\ej) = J p(<p) (ej [ 4>) {4> \ ej) dS n (3) 
2.1 Smeared spectral decomposition 

Usual spectral decomposition of a density operator p = J^Pk l e fe)( e fc| can 
be treated as an atomic measure on CB n whose density is the appropriate 
combinations of the delta functions: 



Pspec((p) = y^Pfc^(0 ~ e fc ) 



For further purposes a 'smeared' version of the spectral decomposition 
is needed, begin with some technical setup. Denote by po the smallest 
eigenvalue of the density matrix p (recall that po > as we restrict ourselves 
to full-range density matrices). Let K be an integer such that K + 1 < 
(1/npo), then the density matrix p is represented as a continuous ensemble 
with positive density: 



^2pk\e k )(e k \ = J n{4>) \4>)(4>\ dS n 



4>eCB n 
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with 



({K+l)n)\ 
K ■ (Kn)\n\ 



1 



(if + l)n 



K e fc I 0)1 



2Kn 



(4) 



Furthermore, the distribution (jlj) tends to the spectral decomposition of p 
as K tends to infinity. See appendix for the proof of formula @. 

2.2 Optimal entropy ensembles 

Let us begin with a single particle case. We need to solve the following 
variational problem. Given a functional Q on L 1 (CS n ) and a density matrix 
p in 7i, find the distribution [i on the set CB n of unit vectors in 7i such that 



r / m 



(5) 



, Q(m) - * ex tr 

We choose the differential entropy of the distribution \x as the functional Q: 



J n(4>) ln/i(0) dS n 



then, according to ©, the variational problem (JSJ) takes the form 



(6) 



— J M^O m M(0) — > extr 

Solving this variational problem by introducing Lagrangian multiples Xij 

we get 



:l+ln /,(</>)) - ^^(e* 



e,} = 



(7) 



Combining the Lagrange multiples into the operator X = X{j |ej)(ej 
we have 



rift 



-(<t>\X\4>)-i 



(8) 
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and the problem reduces to finding \x from the condition 



M0) 10X01 dS n = p (9) 
4>eCB n 

which according to (jHJ) and © and redefining X := — I — X can be written 
as 

(ei\p \ej) = J e^ x W ( e , | 0} (0 1 Bj ) dS n (10) 

<t>eCB n 

It follows from (J7| that the coefficients X^ can be chosen so that X^ = X^. 
That means that the problem of finding the optimal ensemble reduces to 
that of finding the coefficients of a self-adjoint operator, that is, to finding 
n 2 numbers from n 2 equations. 



/ 



2.3 The explicit form for single particle case 

In case of single particle system the equation can be given an explicit 
form. First note that, given a self-adjoint operator X = ^2 x k l e fc)( e fc|> for 
any integrable function f(x) the operator integral 



dS r 



is diagonal in the eigenbasis of X (see the appendix in [Hj for proof). There- 
fore, in order to calculate the integral ()1U|) . we only need its diagonal ele- 
ments. Calculate first the functional K, 4 : 



K{X) = J e <*l*l*>dS n 
for which the following formula holds 

where the eigenvalues of the operator X, 



(11) 



(12) 







1 


1 


1 




. ,x n ) = det 












1 


™.n— 1 
x 2 


„n— 1 
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is the Wandermonde determinant, and the matrix W\ is defined as 



Wi{xi, ...,x n ) 





e Xl 


e X2 . 


e Xn 




1 


1 


1 


det 


XI 


x 2 


x n 




X l 


T n-2 


rj.n-2 



(13) 



We can then explicitly write down the expressions for the functional K 



k=i nj=^ (** *i 



For the operator (|TU)) we have (e^- 1 V/C |ej) 



(14) 



(n-1)! 



+ 



n 1 n 

- xj) n (^fc - n - 

'J i — 1 i — 1 



/ - i A 

fc=l ^ Xfe 

\ Mi / 



(15) 

So, in order to obtain the optimal ensemble for the single particle density 
matrix p = Pj l e j}( e il we solve the system of n equations (e^ | V/C \ej) = 
Py 



3 The bipartite case 

Let p be a robustly separable density matrix in the product space TL®TC ' . 
Then it can be represented (in infinitely many ways) as a continuous ensem- 
ble of pure product states. Carrying out the same reasoning as in section 
12.21 we conclude (see section 13.21 for futher details) that among those contin- 
uous ensembles there exists one having the greatest differential entropy, this 
will be the ensemble we are interested in. Like in section l2~2| formulate the 
variational problem. Let p be a density operator in a tensor product space 
25 = TC Ti! . The task is to find a probability density p.(4>4>') defined on the 
Cartesian product T = CB n x CB n of the unit spheres in 7i, H' , respectively. 



/ 

I Q(m) 



(16) 



extr 
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with 



Q(M) = - J n(<M>')hxn((M>')dS n dS' n (17) 

4>eCB n 

Proceeding exactly in the same way as with single particle, we get the fol- 
lowing representation: 

P= J eWVW) |#'X#1 dS n da' n (18) 

<j><t>'e% 

for some self-adjoint operator X in C. 

3.1 Bipartite separability problem 

From a formal point of view, the equation (|18|) provides a solution of bi- 
partite separability problem. In fact, given a product basis {(ejej'|} in the 
product space 23, (|T%|) is a system of n 4 numeric equations with respect to 
n 4 variables — the matrix elements of the quadratic form X. 

Pii>kk> = J eWIW> (e^ | 00') (00' | e fc e fc ,) dS n dS' n (19) 

If the solution exists, (|19|) provides explicitly an instance of representa- 
tion of p as a mixture of pure product states. 

Although these equations are transcendental. Even in the simple case of 
a single particle in dimension 2, when the operator X is proved to commute 
with p = pi |e 1 )(e 1 1 +p 2 l e 2)( e 2|> we have two variables x±,X2 which we have 
to find from the following system of equations, which are a special case of 
(US): 

' e x i(x 1 -x 2 -l)+e x 2 
(xi-x 2 ) 2 

< 

e x 2(x 2 -x 1 -l)+e x ^ 
- (xi-x 2 )' J 

But the essence of the separability problem is the existence of a solution of 



= pi 

(20) 

= P2 
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3.2 The existence 



In this section I show that the existence of a solution of the equations (|19j) 
for any robustly separable state follows from the concavity of the functional 
Q in (|16|) . Begin with a two-dimensional geometrical analogy. Let S 2 = 
{(xi, X2, x%) I x\ + X2 + £C3 = 1; xi,X2,X3 > 0} be a 2-dimensional simplex 
in R 3 and /C be a functional on the plane P : x\ + X2 + X3 = 1 symmetric 
with respect to x\,X2,xs. If /C is concave, that is, for any a, 6 G P and any 
point g lying between a and b 

K,{q) > IC(a), K{q) > K(b) 

Let K, be bounded in a domain containing the simplex S 2 and let I be a line 
intersecting the interior of S 2 . Then K, has a local maximum in the interior 
of S 2 . 

Now return to the case of bipartite density matrices. Any such matrix 
p can be represented as (an infinite) number of pseudomixtures of pure 
product states (see, e.g. for a discussion). Consider the affine space Pi 
of all normalized signed measures on the Cartesian product T = CB n x CB n 
of unit spheres (it plays the role of the plane P in the example above). Given 
a density operator p, the collection l\ of all pseudomixtures representing p 
is an affine submanifold of Pi. The set of all probability measures on T is 
a simplex S°° in Pi. 

A density matrix p is robustly separable if and only if l\ intersects the 
interior. This because p can be represented as a mixture of products of full- 
range density matrices in Tl, 7i' each of which can be, in turn, represented 
by a smeared spectral decomposition (@J, that means that there exists a 
probability measure representing p which does not vanish anywhere on T 

Therefore the solution of the equations (|19|) always exists for any robustly 
separable p as the functional /C introduced in (|17|) is concave. 

4 Concluding remarks 

First sum up the contents of the paper. Given a density operator p in the 
product of two Hilbert spaces each dimension n, to solve the separability 
problem is to tell if it has a representation by a mixture of pure product 
states. In this paper the most smeared distribution on pure product states 
which yields p rather than their finite mixture is considered. 

The separability problem is reduced to n 4 numerical equations which are 
symbolically written as 



S 



V/C = p 

whose coordinate form in any product basis in 23 = TC (g> H is (fT9")l 



</><£' ex 

These equations are transcendental even for non-product case ()2U(I , but their 
solution always exists for any robustly separable density matrix p. 

That is why one can look for asymptotic methods of finding the solution 
of the separability problem along these lines. This issue will be addressed 
in the next paper on optimal ensembles. 
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A Smeared spectral decomposition 

The formula (j3J) is derived as follows. Let (e| be a unit vector in TC, then 
for any integer m the following formula holds 

/ |(e|^}| 2m |0)(0| dS n = m ) (n ~^ ! (m|e)(e| + I) (21) 
J [m + n)\ 

4>£CB n 

which is a direct consequence of a more general formula obtained in [3] . Let 
m = Kn for some integer K. Then it follows directly from 1)21(1 that 



Kn\e M + l= J { ^ n % \(-\^)\ 2Kn \m dS n 

<t>£CB n 

Dividing this equation by n{K + 1) we obtain the following continuous rep- 
resentation of the projector |e)(e| and the white noise matrix A = I/n: 



e )( e l+ 7? i 7 — = / 77? , ^ /ts m ■ K e l 9)\ mm dS n 



K + l , ' x > K + ln J (K + 1) • (Kn)\n\ 

<t>£CB n 

(22) 

The formula (|22jl is valid for any eigenvector (e&| of p. Form a convex 
combination of the lhs of (|22|) with (yet unknown) coefficients and require 
it to be p: 

When I is replaced by ^ |efe)(e/ c |, the second summand in the lhs of the 
above formula takes the form -^q-j- • ^ • X] l e fc)( e fcl> from which we obtain 

K 1 1 

Ik ■ T7— r + T7— r • - = Pk 
K + 1 a + 1 n 

the we get the expression for ^ 

if + 1 / 1 \ 

qk = ~ir \ Pk - wt^j (23) 

So, we can form the convex combinations of the expressions (|22j) with the 
coefficients qk which yields us p: 
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K V {K + l)nJ J (K + 1) ■ (Kn)\n\ 

k <p£CB n 

Let po be the smallest eigenvalue of p, then, if K + 1 < (1/npo), all the 
coefficients q k in ()23[) are positive. 

B Proof of the formula ( THI ) 

Given a self-adjoint operator X in consider its eigenbasis {e^} and the 
integral 

F n ( Xl ,...,x n ) = J f({<f>\X\<f>))d& n 

CB„ 

In order to calculate it denote (<j) \ e k ) = e ld r k . Then, taking into account 
that the integrand does not depend on phases, the above integral (with 
respect to the normalized measure) reads 



n— 1 \ /n—1 



n-l 

E d<l 



F n (xi, . . . , x n ) = &-jy.(2ir) n [ f[x n + Y,(*k-x n )r 2 k ) [\{r k dr k 

k=i J \k=i 

(24) 



then, introducing the variables t k = ri, we reduce the above integral to 



n—1 \ /n—1 



F n (xi, . . . ,x n ) = (n-l)! / / ix n + ^ (x k - x n ) t k \ dt k 

J \ 1 1 / \ 7 1 



n-l 

E *fe<i 



fc=i / \fc=i 



n-l 



(25) 



For f(z) = e z the above equation reads ^feyN = e Xn / [1 e^ -2 ^ * fe di fc 



n-l fc=l 

fc=i 
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n-2 

i-E h 

k=i 



r n-2 y 1 



E *fc<i 
fc=i 



E *k<i 



= — / Yl e kdt k — - — / Yl e " dt 

n-2 k ~ 1 n-2 fc — 1 

E *fe<i E *k<i 

fc=i fc=i 

and we get the following recurrent relation 



F n (xi,...,x n ) = H — (F n 1 (xi,...,x n -i)-F n 1 (xi,...,x n - 2 ,x n ) 

(26) 

For ii = 2 we have the explicit expression like (|2()|) . for higher n one can 
directly verify that the expression (|14|) satisfies the recurrent relation Q26JI . 
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